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ABSTRACT 

For the initial mass range (140 < M < 260 Mq) stars die in a thermonuclear runaway triggered by 
the pair-production instability. The supernovae they make can be remarkably energetic (up to ~ lO^'^ 
ergs) and synthesize considerable amounts of radioactive isotopes. Here we model the evolution, 
explosion, and observational signatures of representative pair- instability supernovae (PI SNe) spanning 
a range of initial masses and envelope structures. The predicted light curves last for hundreds of 
days and range in luminosity, from very dim to extremely bright L ^ 10'*'* ergs/sec. The most 
massive events are bright enough to be seen at high redshift, but the extended light curve duration (^^ 
1 year) - prolonged by cosmological time-dilation - may make it difficult to detect them as transients. 
An alternative approach may be to search for the brief and luminous outbreak occurring when the 
explosion shock wave reaches the stellar surface. Using a multi-wavelength radiation-hydrodynamics 
code we calculate that, in the rest-frame, the shock breakout transients of PI SNe reach luminosities 
of lO'"^ — 10*^ ergs s^*, peak at wavelengths ~ 30 — 170 A, and last for several hours. We explore the 
detectability of PI SNe emission at high redshift, and discuss how observations of the light curves, 
spectra, and breakout emission can be used to constrain the mass, radius, and metallicity of the 
progenitor. 



1. INTRODUCTION 

Models of metal-free star formation suggest that the 
first stars to form in t he universe were likely quite mas- 
sive, M > 100 Mr^fBromm et al.|[T999 HAbe'l et al.| [2000l: 
iNakamura &: Un icmura 20011 ^ The low metallicity of 
these star may have allowed them to retain much of 
their mass throughout their evolution. The most mas- 
sive objects (M > 260 Mq) are thought to end their 
lives by direct collaps e to a black h ole, with no asso- 
ciated supernova (He ger et al.l |2003D . But stars with 
initial masses between ~ 140 and 260 Mq fall prey to 
the pair-prod u ction i nstability and explo de com pletely 
(IBarkat et al.l 119671: IRakavv et al.l 11967 : Bon d et all 
1984^: 'Umeda fc Nomot ol 120021: IHegenTWoosleyl 120021: 
Waldnian 2008; Moriv a et alTl2010i: iFrver et al.M20l'o[ ). 
The high core temperatures lead to the production of 
e^/e~ pairs, softening the equation of state and leading 
to collapse and the ignition of explosive oxygen burning. 
The subsequent thermonuclear runaway reverses the col- 
lapse and ejects the entire star, leaving no remnant be- 
hind. The explosion physics is fairly well understood and 
can be modeled with fewer uncertainties than for other 
supernova types. 

The predicted explosion energy of pair instability su- 
pernovae (PI SNe) is i mpressive, nearly 10^'^ e rgs for 
the most massive stars (jHeger fc WooslevI 120021) . Ra- 
dioactive ^^Ni can be synthesized in abundance, up to 
40 M0 of it. This is almost 100 times the energy and 
^^Ni yield of a typical Type la supernova (SN la). The 
light curves of the most massive PI SNe are then ex- 
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pected to be very luminous (■^ 10'*^ — 10'*'' ergs/sec) and 
long-lasting 300 days). Deep searches could poten- 
tially detect these events in the early universe, offer- 
ing a means of probing the earliest generation of stars 
(jScannapieco et aLll2005[ ). 

Interest in PI SNe has recently been renewed by the 
discovery, in the more nearby universe, of several su- 
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cases, the high luminosity can be attributed to an in- 
teraction of the supe rnova ejecta with a surrounding cir- 
cumstellar medium (jSmith fc McCravl 120071 ). But other 
events show no clear signatures of interaction, and ap- 
pear to have synthesized large quantities of ^^Ni. The 
speculation is that some of these events represent the 
pair instability explosion of very massive stars, perhaps 
having formed in pockets of relatively lowly enriched 
gas. So far, t he most promising candidate appe ars to 
be SN 2007bi (IGal-Yam et "alll200l IWung et al.l l2010'). 
a Type Ic supernova that was both over-luminous and of 
extended duration. 

In this paper, we model the stellar evolution fi i2.ip and 
explosion fi i2.2p of PI SN models spanning a range of ini- 
tial masses and envelope structures. Using a radiation- 
hydrodynamics code, we then calculate the very lumi- 
nous breakout emission that occurs when the explosion 
shock wave first reaches the surface of the hydrogen en- 
velope f iJ3.1l) . We follow with time-dependent radia- 
tion transport calculations of the broadband light curves 
(i j3.2[) and spectral time-series ( ij3.3p . The models illus- 
trate how the observable properties of PI SNe can be 
used to constrain the mass, radius, and metallicity of 
their progenitors stars. They also allow us to evaluate 
the prospects of discovering these events in upcoming 
observational surveys (Q. 
The detectability of PI SNe was explored previously by 
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TABLE 1 

Pair Instability Supernova Explosion Models 



Name 






A/hc 


b 

Pb 


r„,„kdo''K) 




enersv fB^ 


A/ni 


i?o(10^^cm) 




B150 


150 


150 


66.9 


1.87 


3.51 


68.9 


5.85 


0.00 


4.63 


1.50 


B175 


175 


175 


84.3 


2.31 


3.88 


90.9 


14.6 


0.00 


6.24 


1.82 


B200 


200 


200 


96.9 


2.86 


4.29 


200 


27.8 


1.90 


6.57 


2.13 


B225 


225 


225 


110.1 


3.75 


4.75 


225 


42.5 


8.73 


9.79 


2.46 


B250 


250 


250 


123.5 


5.59 


5.38 


250 


63.2 


23.10 


13.1 


2.79 


R150 


150 


142.9 


72.0 


2.16 


3.70 


142.9 


9.0 


0.07 


162 


1.42 


R175 


175 


163.8 


84.4 


2.66 


4.10 


163.8 


21.3 


0.70 


174 


1.69 


R200 


200 


181.1 


96.7 


3.32 


4.56 


181.1 


33.0 


5.09 


184 


1.76 


R225 


225 


200.3 


103.5 


4.88 


5.15 


200.3 


46.7 


16.5 


333 


2.10 


R250 


250 


236.3 


124.0 


9.45 


6.16 


236.3 


69.2 


37.86 


225 


2.60 


He070 


70 


70.0 


70.0 


2.00 


3.57 


70.0 


8.2 


0.02 






He080 


80 


80.0 


80.0 


2.32 


3.88 


80.0 


17.5 


0.19 






He090 


90 


90.0 


90.0 


2.70 


4.20 


90.0 


28.6 


1.15 






HelOO 


100 


100.0 


100.0 


3.20 


4.53 


100.0 


40.9 


5.00 






HelOOF 


100 


100.0 


100.0 


2.98 


4.44 


100.0 


40.2 


3.64 






HellO 


110 


110.0 


110.0 


4.08 


4.93 


110.0 


55.6 


12.12 






Hcl20 


120 


120.0 


120.0 


5.42 


5.39 


120.0 


70.6 


23.83 






Hol30 


130 


130.0 


130.0 


9.01 


6.17 


130.0 


86.7 


40.32 







this and all other masses in units of Mq 
^ bounce density in 10~® g cm~'^ 
^ presupernova luminosity in 10'*" ergs s~* 

iScannapieco et all ()2005[ ). who used a grey fiux-limited 
diffusion method to model the light curves. In this paper, 
we have generated a new set of more finely resolved mod- 
els which explores the parameter space in a systematic 
way. We have also significantly improved the radiative 
transfer calculations by using a multi- wavelength implicit 
Monte Carlo code which includes detailed line opacities. 
This allows us to generate synthetic spectra and to pre- 
dict the color evol u tion a nd K-correction effects. While 
IScannapieco et all (|2005D focused on the long duration 
light curves, we model here as well the brief and lumi- 
nous transient at shock breakout, and consider whether 
that might be a useful signature for finding PI SNe soon 
after they explode. 

2. EVOLUTION AND EXPLOSION 
2.1. Stellar Evolution Models 

The electron-positron pair-instabilit y supernova mech- 
anism was originally di s cussed by IRakavv fc Shavivl 
(jl967[) and IBarkat et alj (|1967D . and has since been 
explored wit h a number of numeri cal and analytic 
models (see iHeger fc Woosley 120021 and ref erences 
therein). Most recently. Hegcr fc WooslevI ()2002D [here- 
after HW02] explored the explosion of non-rotating bare 
helium cores with a range of masses from 64 to 133 Mq. 
A subset of these models, with masses 70 - 130 Mq (in 
steps of 10 Mq) will be examined here. 

Models of bare helium cores, though computation- 
ally expedient, likely do not fully represent the class of 
PI SNe, as not all stars will have lost their hydrogen en- 
velopes to mass loss or binary mass exchange just prior 
to exploding. We therefore evolved a new set of models 
especially for this study of the breakout transients, light 
curves, and spectraQ They consist of five models each 
of "hydrogenic" (i.e., possessing hydrogen) stars in the 
main sequence mass range 150 to 250 Mq, with initial 
metallicities of and 10~^ times solar. The surface zon- 

^ These models are available to others seeking to carry out sim- 
ilar studies. 



ing of these models was chosen much finer than those 
of HW02 in order to facilitate calculation of the shock 
breakout emission. In addition, one helium core of mass 
100 Mq was calculated with three orders of magnitude 
finer surface zoning (to 10"* Mq) than in IIW02. All 
models were calculated using th e Kepler code a nd th e 
physics discussed in IIW02 and iWooslev et al.l ()2002D . 
One difference with the previous studies is that mass 
loss was included in the hydrogenic stars with non-zero 
metallicity; however, the mass lost was both small and 
very uncertain. 

Properties of all presupernova stars and their explo- 
sions are given in Table 1. The major distinction between 
the zero metallicity and 10^** solar metallicity presuper- 
nova stars was that the former died as compact blue su- 
pergiants (BSG) while the latter were red supergiants 
(RSG) with radii 10 to 50 times larger. To some extent, 
this difference also relies on a particular choice of un- 
certain parameters - in particular primary nitrogen pro- 
duction and mixing - so not too much weight should be 
placed on the metallicity of the model. The production 
of primary nitrogen was, by design, minimal in all mod- 
els. In the zero metallicity series, especially the higher 
mass ones, this required reducing semi-convection by a 
factor of 10 compared with its usual setting in Kepler 
and turning off overshoot mixing. Unless this was done, 
only the 150 Mq model ended up as a blue supergiant 
while the other four were red owing to primary nitrogen 
production. 

For the 150, 175, and 200 Mq stars with 10""* solar 
metallicity, both semi-convection and overshoot mixing 
had their nominal settings, but red supergiants resulted 
without primary nitrogen production. The 225 and 250 
Mq stars with 10"'' solar metallicity used the reduced 
mixing - again to avoid large ^^N production - but also 
ended up as red supergiants. Since even a moderate 
amount of rotation would lead to some mixing between 
the helium core and hydrogen envelope in the zero metal 
stars, and since the nominal semi-convection and over- 
shoot parameters also lead to primary nitrogen produc- 
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Fig. 1. — Ejecta density profiles of PI SN models 10 days after explosion, when the expansion is nearly homologous [v oc r). 




Fig. 2. — Ejecta abundance structures of three of the PI SN models in the homologous expansion phase. An artificial smearing of the 
compositional interfaces has been applied to mimic multi-dimensional mixing processes. 



tion, it is likely that a fraction, perhaps all, of the zero 
metallicity pair instability stars also die as red super- 
giants. This would imply that blue supergiants are a 
rare population of supernova progenitors for such mas- 
sive stars even at ultra-low metallicity. Our goal here 
though was to prepare a set of blue and red supergiant 
progenitors with a range of masses to explore how the ob- 
servable properties of the explosion may depend on the 
structure of the hydrogen envelope. 

2.2. Explosion 

The pair instability is triggered in helium cores above 
about 40 M0 once the temperature in the stellar core 
exceeds ~ 10^ K, i.e., after helium burning and during 
carbon ignition. The creation of e"'"/e~ pairs in the core 
then softens the equation of state to below 7 < 4/3 lead- 
ing to instability and collapse. For helium cores below 
about 135 Mq the collapse is eventually reversed by ex- 



plosive nuclear burning, first of oxygen and, for more 
massive cores, silicon. Non-rotating helium cores more 
massive than 133 Mq experience a photo-disintegration 
instability after silicon burning and collapse directly to 
a black hole (HW02). Between 40 and 60 Mq, the pair 
instability in bare helium cores leads to multiple violent 
mass eje ctions, but does not d isrupt the entire star on the 
first try (jWooslev et al.ll2007[ ). Between 60 and 133 Mq, 
helium stars collapse to increasing central temperature 
and density, explode with greater violence, and produce 
more heavier elements, especially ^^Ni (see Table [T]). 

Most of the helium core explosions studied here were 
taken from the survey by HW02, but the 100 Mq model 
was recalculated with finer surface zoning. A possible 
point of confusion is whether shocks ever form in these 
sorts of helium stars or whether, given their small ra- 
dius, the surface remains in sonic communication with 
the center throughout the collapse and initial expansion. 
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Fig. 3. — Calculated bolometric light curves of the shock breakout transient from pair-instability explosions of massive red-supergiant 
(left) and blue-supergiant (right) stars. The duration of the burst, determined by the light crossing time, is significantly longer for the 
more extended red supergiant models. 



We find that for tlie 100 Mq model the outer layers do 
not participate in the collapse and, at about 2 seconds 
after maximum compression, a very strong shock forms 
initially about 2 Mq beneath the surface at a radius of 
8 X 10^ cm. Due to the large explosion energy and ac- 
celeration in the steep density gradient, material in the 
shock reached a speed of about one-third the speed of 
light before erupting through the photosphere. 

The explosions of the hydrogenic stars had character- 
istics, including kinetic energy and nucleosynthesis, ap- 
proximately set by the mass of their helium cores (Table 
[T]). There were, however, some important differences. 
For the helium stars the mass was set at the beginning 
of the calculation and held fixed, while in the hydrogenic 
stars the helium core grew significantly after central hy- 
drogen depletion owing to hydrogen shell burning. In 
some cases it also shrank due to convective dredge up. 
At death the helium core thus had a different composition 
and entropy than a corresponding helium star evolved at 
constant mass. 

An even more significant difference of the hydrogenic 
stars is that their exploding helium cores encounter the 
lower density hydrogen envelope and interact with it 
hydrodynamically. This has several important conse- 
quences. First, the expansion of the helium core is 
slowed, which produces Rayleigh- Taylor instabilities and 
mixing. Up until this point, the explosion had been de- 
termined by simple physics and (assuming no rotation) 
well represented by a spherical one-dimensional calcula- 
tion. While the mixing_can be calculated in a multi- 



2009bD. it has not 



dimensional code (e.g., Jogg erst et al.l ! 
been y et and is parametrized here as in[ 
(j2009( ). Except for two models that experienced signif- 
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icant fallback (see below), the mixing has no effect on 
the nucleosynthesis or energetics and does not affect the 
breakout emission. However, the late time spectra may 
be somewhat sensitive to mixing. Figures [T] and [2] plot 
the final density and mixed compositional structures for 
a few representative models. 

In some cases, the deceleration of the helium core by 
the envelope is so severe that the star doesn't completely 
explode. This is particular ly true for the b lue supergiant 
m odels. As pointed ou t bv lChevalied (I1989D and ex plored 
bv lZhang et all (|2008D and lJoggerst et al.l (|2009al) . brak- 
ing and fall back occur to a greater extent in more com- 
pact stars. The helium core encounters its own mass 
earlier and the reverse shock returns to the center when 
the density is still high. In our case, this resulted in mod- 
els B150 and B175 failing to eject anything other than 
their hydrogen envelope. Of course the story does not 
end there ~ these bound helium cores oscillate awhile, 
settle down and evolve again. Given their residual he- 
lium masses, models B150 and B175 will probably be- 
come pair instability supernovae again and disrupt en- 
tirely on the second try. They thus represent an exten- 
sion to higher masses of the pulsational p air instability 
supernovae studied bv fWooslev et all (^7). Depending 
upon the timing, the second mass ejection may produce 
an extremely bright supernova as the very energetic sec- 
ond supernova plows into the ejecta of the first. This 
future evolution is beyond the scope of the present pa- 
per, but the evolution and explosion of blue supergiants 
with ZAMS masses of 100 - 200 Mq is clearly an area 
worth further study. 

3. OBSERVABLE PROPERTIES 
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Fig. 4. — Synthetic spectra of the shock breakout transients averaged over the burst peak. The time-averaged spectra are the convolution 
of several blackbodies at different temperatures, and so are broader than a single blackbody function peaking at the same wavelength 
(dashed lines). Due to the low thermalization opacity, the color temperature is greater than the photospheric temperature by a factor of 
2-3. 



3.1. Shock Breakout 

In the hydrogenic models, the expansion of the ex- 
ploded helium core into the hydrogen envelope drives 
a radiation dominated shock. When this shock ap- 
proaches the surface of the star, the postshock ra- 
diation can escape in a luminous x-ray/UV burst 
jKlein & Chevalier 119781 : lEnsman fc Burrows! 119921 : 
iMatzner fc McK ee 1999*). This event occurs at a distance 
AR from the stellar surface such that the diffusion time 
from the shock front {td ^ tAR/c) is comparable to the 
dynamical time for the shock to travel the same distance 
{te ^ AR/vg). This implies an optical depth r w c/vg 
at breakout. The shock velocity Vs scales with (E/MY^^ 
and for pair SNe is of order 3 x 10'^ — lO"* km s~^, sim- 
ilar to that of ordinary core collapse supernovae. Shock 
breakout thus occurs at optical depths of Tf, « 30 — 100. 

The surface layers of PI SNe progenitors are difficult to 
model; a proper representation of the atmospheric struc- 
ture would require a detailed radiation transport within 
the stellar evolution code and possibly the inclusion of 
3-D effects. The limitations of the 1-D stellar evolution 
models thus introduce some uncertainty into the break- 
out predictions. While the current models were much 
more finely zoned than previous calculations, they still 
did not fully resolve the optically thin layers of the star. 
They did, however, determine the slope of the steep den- 
sity profile at the surface, out to an optical depth of a 
few. To extend the profile into the optically thin region, 
we fit a power law to the outermost zones and linearly 
interpolated and extrapolated the density structure. The 
re-zoned model had 100 zones within the region t < Tt, 
and extended to a minimum r of 10~^. 

We followed the supernova explosion using the Kepler 
code until the shock front began to approach the stel- 
lar surface (r w 10^). The model structure was then 



mapp ed into a modified version of the SEDONA transport 
code (|Kasen et al.ll2006l) which coupled multi- wavelength 
radiation transport to a staggered mesh 1-D spherical 
Lagrangian hydrodynamics solver. A standard artifi- 
cial viscosity prescription was included to damp oscil- 
lations behind the shock front. We followed the ra- 
diation transport using im plicit Monte Carlo methods 
(jFleck fc Cummin"gill971[ ) coupled to the hydrodynam- 
ics in an operator split way. The transport was treated 
in a mixed-frame formalism, in which photon packets 
were propagated in the inertial frame, but the opacities 
and emissivities were computed in the comoving frame, 
and the proper Lorentz transformations were applied to 
move between the frames. The radiation energy and mo- 
mentum deposition terms in the hydro equations were 
estimated by tallying the energy and direction of packets 
moving through each zone. Further details of the nu- 
merical methods in the contex t of shock breakout will b e 
given in a separate publication iKasen fc WooslevI ()201lD . 

The opacity and emissivity of the models were dis- 
cretized into 25000 wavelength bins covering the range 
A = 0.1 — 25000 A, while the output spectra were binned 
to a coarser resolution to improve the photon statis- 
tics. The dominant opacity in the supernova shock is 
electron scattering, while free-free is the most impor- 
tant absorptive opacity. Because the current models 
had zero or very low metallicity, we ignored bound-free 
and line opacity. At wavelengths near the blackbody 
peak, the ratio of free-free opacity to total opacity is 
typically small Ea ^ 10""*. In the postshock region, 
Compton up-scattering may then become an important 
mean s of e nergy exchange between matter and radia- 
tion d^^er 1976). On average, the fractional change 
of a photon's energy in a single Compton scattering is 
Cc = AkT/nieC^ , which gives an effective ec — 6 x 10~^ 
at T = 10^ K. Comptonization will significantly alter 
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TABLE 2 
Shock Breakout Transients 



name Lp^^k (ergs s"!) At,,'' T,ff(10^) T,oi(105)'= Ap'^ Ep'^ Etot" E^lJ 



R250 9.6xl0*-^ 5858 1.3 3.5 82.5 0.15 5.9x10-*^ 2.2x10''^ 

R200 2.6x10*^ 6422 1.0 2.2 130.5 0.09 1.8x10-*^ l.exlO'"' 

R150 1.2x10*5 7051 0.9 1.7 168.5 0.07 9.3x10-*** 1.3x10*'' 

B250 1.4x10*5 966 3.3 6.3 45.5 0.27 1.3x10** 6.5x10** 

B200 6.7x10** 714 3.9 7.7 37.5 0.33 5.8x10*^ 1.9x10** 



^ duration of burst, full-width half maximum, in seconds 
^ color temperature 

spectral wavelength peak, in A 

spectral energy peak, in keV 
" total energy emitted in burst, in ergs 

' total energy emitted in burst at wavelengths greater than 

the radiation spectrum after scattering through an op- 

— 1/2 

tical depth Tc ^ €c which, for higher temperatures, 
may be less than the thermahzation depth to free-free 

— 1/2 

opacity Ta — ea ■ In our calculations, we therefore ap- 
proximated the Compton effects by taking a fraction Ec 
of the Thomson opacity to be thermalizing. This grey 
absorptive component was added to the wavelength de- 
pendent free-free opacity. Clearly a direct treatment of 
inverse Compton scattering is needed to more accurately 
predict the spectrum of the breakout burst; nevertheless, 
our treatment of the full non-grey radiation transport of- 
fers some advance over previous numerical calculations. 

Once a shock front reaches the steep outer layers of 
the star, it accelerates down the steep density gradient. 
For model R250, the shock velocity reaches w 1.5 x 
10^ km s^^ around breakout. The temperature of the 
post-shocked gas is determined by the jump conditions 
= pqv'^ / + 1) , where po is the density of the pre- 
shocked material and 7 = 4/3 is the adiabatic index for 
a radiation dominated gas. Near the surface, the RSG 
models had densities of po ^ 10~^^ g cm~'^ which gave 
typical postshock temperatures of Tg « 5 x 10^ K. The 
BSC models, being more compact, had higher densities 
at the surface and hotter temperatures, Tg « 10^ K. 

Figure [3] plots the calculated bolometric light curves of 
the breakout transients. The key observable properties 
are listed in Table [2J These calculations included the 
proper light travel time effects, which primarily deter- 
mined the duration of the observed burst. At w Rq/c. 
The RSG models lasted 1 to 2 hours, and reached peak 
luminosities of Lpeak = 10"*^ - lO'*^ ergs s'^. The BSC 
models, being more compact, had briefer transients last- 
ing 0.1 hours, with somewhat lower peak luminosities, 
-^poak ~ ergs/sec. These values can be compared to 
the breakout in an ordinary Type IIP supernova explo- 
sion: Lpoak = 5 X 10"'"^ ergs s~\ AT — 0.5 hours. As 
it turns out, the shock velocities and post-shock energy 
densities of PI SNe are comparable to that of Type IIP 
supernovae, but the breakout emission can be signifi- 
cantly brighter due to the larger stellar radii. 

Figure 0] shows synthetic spectra of the bursts, aver- 
aged over the peak of the breakout light curve. The 
RSG spectra peak at wavelengths Ap = 80 — 170 A 
(~ 0.07 — 0.15 keV), while the BSC spectra peak near 
37 — 45 A 0.3 keV). These peak wavelengths are com- 
parable to that of ordinary Type IIP breakout (Ap « 



Lyman-a, in ergs 

100 A). The spectra are reasonably approximated by a 
blackbody, but with excess emission at both high and 
low frequencies. This is because the time-averaged spec- 
tra represent a convolution of several blackbodies at dif- 
ferent temperatures. 

Because the opacity in the atmosphere is strongly 
scattering dominated, radiation is thermalized at opti- 

— 1/2 

cal depths t w ea well below the photosphere. The 
emergent spectrum is thus characterized by the higher 
temperature of these deeper layers. Defining the color 
temperature Tcoi of the mean spectrum as the temper- 
ature of a blackbody peaking at the same wavelength, 
we find that the observed Tcoi is a typically a factor 
2-3 times higher than the effective temperature at the 
photosphere (Table 2). A similar effect was noted by 
fEnsman fc Burrowj ()1992[) for m odels of SN 1987A . In 
addition, recent analytic work by iKatz et all ()2010[) has 
shown that when the shock velocity is relatively high, 
> 0.1c (and when the thermahzation by Compton scat- 
tering is treated properly) non-equilibrium effects may 
lead to a significantly harder non-thermal component in 
the spectrum. In the BSC models, the shock velocities do 
in fact reach ~ 0.1c and such a non-thermal component 
could modify the predicted spectra shown here. The the 
RSG models, on the other hand, have lower shock veloc- 
ities and are likely not strongly affected. 

Following shock breakout, radiation continues to dif- 
fuse out of the expanding, cooling ejecta. This leads to 
a longer lasting emission that may be visible at optical 
wavelengths for several weeks. Such an initial thermal 
component to the light curves is discussed in the fol- 
lowing section. The extreme luminosity of the breakout 
bursts themselves, and their short durations, make them 
appealing transient to search for in high-redshift surveys. 
We consider their detectability in Section |4l 

3.2. Light Curves 

Following shock breakout, the luminosity of PI SNe can 
be powered by three different sources: (i) The diffusion 
of thermal energy deposited by the shock; (ii) the energy 
from the radioactive decay of synthesized ^^Ni; or (iii) 
the interaction of the ejecta with a dense surrounding 
medium. The internal energy suffers adiabatic losses on 
the expansion timescale tex = Ro/v, so source (i) will be 
most significant for stars with large initial radius Rq. In- 
teraction has not been included in the models discussed 
here, but as already mentioned, can be very significant 
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Fig. 5. — Bolometric light curves of the full set of models representing the explosion of blue supergiant stars (left panel), red supergiant 
stars (middle panel) and bare helium cores (right panel). The more massive stars have brighter and longer lasting light curves. 

equilibrium, typically a reasonable approximation for su- 
pernovae in the photospheric phase. 

The resulting model light curves (Figure [5]) span a 
wide range of luminosities and durations. The more mas- 
sive explosions are bright for over 300 days, with lumi- 
nosities exceeding lO'*'* ergs/sec. The long duration re- 
flects the timescale for photons to diffuse through the 
optically thick supernova ejecta. In a homologously ex- 
panding medium, the effective diffusion time scales as 

td - Ki/^Mf/^£:- i/4, where K is the effective mean opac- 
ity (|Arnettiri980D . Model He 130, for example, is nearly 
100 times as massive and energetic as a SN la, and thus 
has a light curve '^10 times as broad. The slow release 
of energy moderates the emergent luminosity, so that 
despite making over 60 times as much ^^Ni as a typical 
SNe la, model IIel30 is only ^ 10 times brighter than 
one at peak. 

The morphology of the model light curves depends on 
the envelope of the progenitor star. The RSG mod- 
els resemble Type II plateau supernovae (SNe IIP) light 
curves, though with luminosities and durations 3 times 
greater. The initial luminosity is powered by the ther- 
mal energy in the extended hydrogen envelope. Over 
time, as the ejecta cool, a hydrogen recombination front 
propagates inward in mass coordinates, which increases 
the ejecta transparency due to the elimination of elec- 
tron scattering opacity. This effect regulates the re- 
lease of the thermal energv ([Grasberg &: Nadezhinlll976l : 
IPopovl 119931 : IKasen fc Wooslevl l2009f ). In the inner re- 
gions, heating from radioactive decay delays recombi- 
nation and causes the light curve to rise to a peak at 
around 200 — 300 days. Eventually, the transparency 
wave reaches the base of the hydrogen envelope, af- 
ter which recombination proceeds much more rapidly 
through the heavier element core. After a rapid release 
of the remaining internal energy, the light curve drops 
off sharply and follows directly the radioactive energy 
deposition rate. 
The light curves of the BSG models have a dimmer 



Fig. 6.— Synthetic R-band light curves (at z = 0) of bright PI SN 
models - R250 (dashed-dot), B250 (solid), and Hel30 (dashed) - 
compared to observations of a normal Type la supernova SN 2001el 
(red triangles, Krisciunas et al. 2003) a normal Type IIP super- 
nova SN 1999em (blue squares, Leonard et al..ii2002^ and the over- 
lumin ous core-collapse event SN 2006gy (green circles. fSmith et al.l 
[2009 ). 

for stars which undergo pulsations before exploding com- 
pletely. 

We calculated light c urves of the e xplosion models us- 
ing the SEDDNA code (|Kasen et al.l [2006). The initial 
density, composition, and temperature structures were 
taken from the Kepler calculations extended into the 
nearly homologous expansion phase (~ 10 days after ex- 
plosion). The energy deposition from ^^Ni decay was fol- 
lowed using a multi-wavelength transport scheme treat- 
ing the emission, propagation, and absorption of gamma- 
rays. For the optical radiation transport, the opacities 
used included electron-scattering, bound-free, free-free, 
and the aggregate effect of millions of Doppler broad- 
ened line t ransitions treated in th e expansion opacity 
formalism ([Eastman fc Pintolll993[ ). Atomic level pop- 
ulations were calculated assuming local thermodynamic 
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Fig. 7. — Broadband light curves of PI SN models. The light curves displayed in each panel are, from bottom to top, UBVRIJHJK. The 
y-axis corresponds to the B-band light curve, while the other light curves are offset by 1 magnitude from the one above or below. 



initial thermal component, due to the relatively smaller 
radii of the progenitors. In this respect they resemble 
scaled up versions of SN 1987A. In models B200, B225, 
and B250, the light curves rise to a bright ^^Ni pow- 
ered peak at ~ 300 days after explosion. The less mas- 
sive events (B150 and B175), which failed to eject any 
^^Ni, only show the initial thermal light curve component 
lasting ^ 150 days, and more closely resemble typical 
SNe IIP. 

The helium core PI SNe models, being the most com- 
pact progenitors, lack a conspicuous thermal light curve 
component altogether. The lower total mass leads to rel- 
atively briefer light curves, peaking around 150 days after 



explosion. Model Hel30 40 Mq of ^^Ni) reaches an 
exceptional peak brightness of 2 x 10^'* ergs s~^ at around 
180 days. Model HeTO, on the other hand, proves that 
despite being massive and energetic, not all PI SNe are 
bright. This explosion produced only 0.02 Mq of ^^Ni 
and the light curve reached only 3 x 10^^ ergs/sec, rather 
sub-luminous for a supernova. 

Figure [6] compares the synthetic R-band light curves 
of representative PI SN models with observatio ns of a 
typic al Type la supernova (SN 2001el, Krisciuna s et al.l 
20031) ■ a tvp i cal T ype IIP supernova (SN 1999em, 
Leonard et al]|2002D . and SN 2006gy, one of the most 



the luminous core-collapse SNe discovered (jSmith et al.l 



Pair Instability Supernovae 



9 




Fig. 8. — Time evolution of the velocity measured at the electron scattering photosphere for several PI SN models. In the BSG and helium 
core models, the photosphere initially recedes quickly as the outer layers recombine and become transparent. The rise in photospheric 
velocity seen in some models reflects the outward diffusion of radioactive energy which may heat and reionize the external layers. 
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Fig. 9. — Spectra of select models at the peak of the light curve. 
The important line absorption features are marked. 

I2006f ). The later has been suggested to be a pair instabil- 
ity supernovae, however the predicted model light curve 
durations are seen to be too long even for this event by a 
factor of several. Another possibility is that SN 2006gy 
was a n example of a puls ational-pair instability super- 
nova (jWooslev et al.l[200l . 

Figure [7] shows the multi-color optical and near- 
infrared light curves for the models. At bluer wave- 
lengths, the light curves generally peak earlier and de- 
cline more rapidly after peak. This reflects the progres- 
sive shift of the spectral energy distribution to the red 
over time. This shift is due not only the decrease in effec- 
tive temperature, but also to the increase in line opacity 
from iron group elements of lower ionization stages (Fe II 
and Co II), which blankets the bluer wavelengths. 

The brightest helium core models such as HelSO dis- 



play a pronounced secondary maximum in the infrared 
light curves. This feature is similar to what is observed 
in Type la supernov ae and the p hysical explanation is 
essentially the same ()Kasenll2006D . When the tempera- 
ture in the ejecta drops below ~ 7000 K, doubly ionized 
iron group elements begin to recombine. As the infrared 
line emissivity is much greater for singly ionized species 
(Fe II and Co II) the flux can be more efficiently redis- 
tributed from bluer to redder wavelengths, leading to an 
increase in the infrared luminosity. The secondary max- 
imum is therefore more prominent in models with large 
abundances of iron group elements. Detection of a sec- 
ondary maximum in an observed supernova would pro- 
vide strong evidence that the explosion did indeed syn- 
thesize large amounts of ^^Ni. The lack of a secondary 
maximum does not necessarily rule out the presence of 
substantial ^^Ni, as strong radial mixing of the nickel can 
sometimes smear out the two bumps ((Kascn.2006.) . 

3.3. Spectra 

The spectra of the PI SN models (Figures[9]and [TU|) re- 
semble those of ordinary SNe, with P-Cygni line profiles 
superimposed on a pseudo-blackbody continuum. Be- 
cause of the low abundance of metals in unburned ejecta, 
many of the familiar line features are weak or missing in 
the early time spectra. For example, at maximum light 
models R250 and B250 show only features from the hy- 
drogen Balmer lines and calcium in their spectrum. How- 
ever, at later times, when the photosphere has receded 
into layers of burned material, other line features appear. 

The maximum light spectrum of the bare helium core 
model HelOO is dominated by lines from freshly syn- 
thesized intermediate mass elements, and resembles a 
Type Ic SN with lines due to Mg I, Mg II, Si II, Ca II 
and O I. After peak, the spectrum shows more features 
from the iron group elements in the ^^Ni-rich core. He- 
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Fig. 10. — Spectral evolution of model R250 (left panel) and 
model HelOO (right panel). Over time, additional line features 
appears as the photosphere recedes into layers of burned material. 
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Fig. 11. — The day 100 spectrum of model R150 shown for dif- 
ferent values of metallicity in the hydrogen envelope, from top to 
bottom Z = 10~^, 10~*, 10~^, and 1 times solar. The metallicity 
can be probed observationally by examining metal line absorption 
features or the flux in the ultra-violet. 

Hum lines are not present at any epoch, as the envelope 
temperatures are too low to thermally excite the lower 
atomic levels of the optical transitions. However, if ^^Ni 
is mixed out into the helium rich layers, non-thermal 
excitation by radioactive decay products c ould generate 
significant helium line opacity (|Lucvl[l99ll ). 

Although PI SNe are highly energetic events, their 
large ejected masses imply only moderate characteristic 
velocities: v ~ (2£;K/Mej)i/2 - 5000 km s"\ about half 
that typical of SNe la and several times less then the 
broad lined Type I c SNe that have been associated with 
gamma-rav bursts (|Galama et al.|[l998[ ). Figure [5] show 
the time evolution of the velocity measured at the elec- 
tron scattering photosphere. In the BSG and helium core 
models, the photosphere initially resides in the outer. 
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Fig. 12. — The spectrum, observed near peak, of SN 2007bi (red 
line Gal- Yam et al. 200(|) compared to the synthetic spectrum of 
our model HelOO at 100 days after explosion (black line). Line iden- 
tifications for the model are marked. Most of the major observed 
line features are reproduced by the model, except for the forbidden 
Ca II line emission at 7300 A which is the result of non-equilibrium 
effects not included in the calculations. 

high velocity layers of ejecta, but recedes quickly as these 
layers recombine and become transparent. Eventually, 
the photosphere settles in the inner regions of burned 
material, where radioactive energy deposition maintains 
the ionization state. Interestingly, the photospheric ve- 
locity in model B250 in fact increases for some time pe- 
riod as the radioactive energy diffuses outward and reion- 
izes the hydrogen envelope. In contrast, the decline in 
photospheric velocity in the RSG models is gradual until 
the recombination front reaches the base of the hydrogen 
envelope, after which it drops off sharply. 

Certain features in the spectra of PI SNe could, if ob- 
served, offer direct confirmation that the progenitor star 
was of low metallicity. The metallicity has two principle 
spectroscopic effects (Figure fTTj) . First, for metallicities 
Z > 10^^, some absorption lines from intermediate mass 
elements are noticeable, for example those of Ca II H&K 
(near 3800 A) and the IR triplet (near 8300 A). Second, 
higher metallicity leads to a significant reduction in the 
ultraviolet (UV) flux, due to the iron group line blan- 
keting at bluer wavelengths. Similar metallicity effects 
have been noted in the spectral modeling of Type IIP 
SNc (Baron et al. 2003; Dcssart & HiUicr 2005) Spectro- 
scopic or rest frame UV observations of PI SNe may, 
together with modeling, constrain the metallicity of the 
progenitor star. This assumes, of course, that the super- 
nova has been observed early enough that we are seeing 
layers of ejecta unaffected by the explosive nucleosynthe- 
sis. 

To draw some comparison to observations. Figure [T^] 
shows the near maximum light spectrum of one model, 
HelOO, compared to the most promising candidate 
to date for an obser ved pair supernova, SN 2007bi. 
iGal-Yam et alJ (|2009l ) has previously shown that the 
light curve of this supernova is a good match to model 
HelOO. On the whole, the correspondence of the spectra 
is also rather good, especially considering that the model 
has not been tuned to match the data. Most of the major 
line features are reproduced, in particular the unusually 
prominent Mg I and Mg II lines. The ejecta velocities (as 
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Fig. 13. — Observer frame K-band light curve of model R250 as a 
function of redshift. The effects of cosmological redshift, dimming, 
and time-dilation have all been included. For 2 > 7, one observes 
in the rest-frame UV, and the initial thermal component of the 
light curve is brighter than the later radioactively powered peak. 



5.0 




years since explosion 

Fig. 14. — Observer frame R-band light curve of model R250 as a 
function of redshift. The effects of cosmological redshift, dimming, 
and time-dilation have all been included. 

evidenced by the blueshift of the absorption features) are 
also roughly in agreement. The model fails to reproduce 
the forbidden Ca II line emission at 7300 A, but this fea- 
ture results from non-equilibrium effects which are not 
included in our calculations. 

4. DETECTABILITY 

Some PI SNc should be observable out to large dis- 
tances, but the effects of cosmological redshift and time- 
dilation effects will significantly affect the shape and lu- 
minosity of the observed light curve. Figures [T3] and [H] 
demonstrate these effects for model R250 in the K- and 
R-bands, respectively. The radioactively powered optical 
emission, which lasts for hundreds of days, emits very lit- 
tle fiux at rest wavelengths A < 2000 A. Thus, to follow 
objects at z > 5 it is best to observe at infrared wave- 
lengths. Alternatively, one could search for the brief and 
very blue emission from shock breakout. 

In Figures [15] and [16] we plot the observer frame R- and 



K-band peak magnitudes of several PI SN light curves 
(starting t > 10 days after explosion in the rest-fame) as 
a function of redshift. Future R-band surveys, such as 
with the Large Synoptic Survey Telescope, will routinely 
reach limiting magnitudes of Mj^ 24.5, which would 
detect or place constraints on the brightest PI SNe out 
to redshift of z ~ 2. Deeper searches (M^ ~ 28) would 
reach redshifts z ~ 5, while a similar K-band search 
would probe redshifts of z ^ 10 and beyond. Space 
based observations at wavelengths longer than K-band 
could potentially see to even greater distances. The long 
duration of the PI SN light curves - greatly prolonged 
by the (1 -t- z) cosmological time dilution factor - poses 
a challenge for detecting them as transients. At z ^ 7 
the light curve of a PI SN can last 1 to 5 years in the 
observers frame (Figures [T3] and [Ti)) . Deep exposures 
over long time baselines would be needed to discover and 
follow such an event. 

An alternative way of discovering PI SNe would be to 
search for the short-lived shock breakout transient. At 
z = 1, the observed burst would last a few hours, and 
at z = 10 about a day. The spectrum, which peaks in 
the rest-frame at A ~ 80 — 170 A, would be seen on 
the Rayleigh- Jeans tail. Unfortunately, neutral hydro- 
gen along the line of sight likely absorbs all radiation 
shortwards of the Lyman alpha line at 1215 A, which 
is the bulk of the radiation emitted in shock breakout. 
However, the radiation at longer wavelengths can be sig- 
nificant and may still be visible. 

Figure [T7I plots the observed AB magnitudes (in differ- 
ent wavelength bands) of the RSG models at the peak of 
the bolometric breakout light curve. At redshifts z > 3 
the curves flatten out as flux is redshifted into the ob- 
served band and compensates for the greater distance. 
An all sky optical survey like LSST which reached mag- 
nitude R « 24.5 per image could potentially catch shock 
breakout to redshifts of z w 1. In principle, deeper opti- 
cal imaging in select fields could detect breakout to much 
higher redshift, but in practice, the Lyman limit restricts 
the optical detectability to z ~ 6. 

Near-infrared surveys could avoid the problem of hy- 
drogen absorption. A facility like JWST with < 10 nJ 
{Mab ~ 29) sensitivity in the wavelength region 2 — 5^ 
has the capability of detecting shock breakout out at red- 
shifts of z « 10 — 20. However, the small field of view of 
JWST, combined with the expected low rates of PI SNe, 
suggest that it will be unlikely to discover a PI SNe in 
this manner. A wide-field imager with infrared capa- 
bilities (e.g., WFIRST or EUCLID) may fare better for 
discovering shock breakout at the highest redshifts. 

The rates of PI SNe, in eit her the nearby or d i stance 
universe, are very uncertain . iScannapieco et aT] ()2005f ) 
(see also IScannapieco et all (|2003l )) provided one esti- 
mate of the rate as of one PI SN per 1000 solar masses 
of metal free stars, with the metal free stars themselves 
forming, at all redshifts, at only 1% of the total star for- 
mation rate. Observationally, the detection of at least 
on promising candidate pair SN (SN 2007bi) suggests a 
rate i n the local universe of 5, 10~^ of the core collapse 
rate (jQuimbv et all [20091 iGal- Yam et al.ll2009f) . More 
detailed investigation of the PI SN rate will be consid- 
ered elsewhere (Pan et al., in prep) but in any case it 
should be kept in mind that PI SNe are rare compared 
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Fig. 15. — Detectability of pair instability supernovae light 
curves. The figure plots the observer-frame peak K-band mag- 
nitude of several models as a function of redshift. 
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Fig. 16. — Same as Figure [TSl but for the R-band. 

to other supernovae, and will therefore be difficult to dis- 
cover in bulk. 

5. SUMMARY AND CONCLUSIONS 

We have surveyed the spectra, multi-color light curves, 
and shock breakout transients of pair instability super- 
nova models throughout the mass range in which the 
star suffers complete disruption. Three varieties of pro- 
genitor stars were explored: i) helium cores that may 
have resulted when very massive stars either lost their 
envelopes to some sort of non-radiative mass loss (as in 
Eta Carina) or to a binary companion; ii) blue super- 
giants still retaining most of their hydrogen envelope; and 
iii) red supergiants with enormous radii. The frequency 
and distribution with mass and metallicity of ii) and iii) 
depends on uncertain parameters of semiconvection and 
overshoot; in general it seems that the red supergiants 
are a more common outcome of the stellar evolution than 
blue supergiants. 

Though PI SNe are often considered to be categor- 
ically bright events, in fact it is only the most mas- 
sive stars (M > 200 Mq) that are extraordinary in 



terms of luminosity. The lower mass objects have peak 
brightnesses comparable to normal core collapse events 
L = 10*^ — lO"*^ ergs, though the light curves last for 
much longer. For a sensible initial mass function, it is 
likely that these dimmer events are actually the more 
representative PI SNe. An extended light curve, rather 
than an extreme luminosity, may therefore be the most 
relevant signature of the pair instability events. 

The models illustrate how the different classes of pro- 
genitor stars can be distinguished observationally by the 
light curve morphology. The red supergiant explosions 
have light curves with long duration plateaus, similar to 
ordinary Type IIP supernovae but with luminosities and 
durations ^ 3 times greater. The light curves of the blue- 
supergiant explosions more closely resemble SN 1987A, 
with a brief initial thermal component and a late, promi- 
nent ^^Ni powered peak. The bare helium core model 
light curves resemble very long lasting Type Ib/Ic super- 
novae. 

The exploding cores of the more compact blue su- 
pergiant progenitors are braked extensively by hydrody- 
namic interaction with their envelopes. For the lighter 
stars considered (models B150 and B175) the system 
failed to completely explode on the first try, and only 
the outer envelope of the star was ejected. The class of 
pulsational pair instability ma.y thus have a larger range 
than stated by iWooslev et al.l ()2007[ ). A secondary ex- 
plosion of the helium core, and its subsequent collision 
with the ejected hydrogen envelope, could then lead to 
an extremely luminous supernova event. 

Despite the large explosion energy of PI SNe, the ex- 
pansion velocities, as measured from the Doppler shift of 
spectral lines, are rather low v « 5000 — 10000 km s^^. 
The spectra of helium core explosions resemble Type Ic 
supernovae, and in particular are distinguished by promi- 
nent lines of Mg I and Mg II. The spectra of hydrogenic 
models appear in most respects like ordinary Type II su- 
pernovae, however the absence of certain metal line fea- 
tures (e.g., the Ca II IR triplet) and the bright ultraviolet 
emission (due to the reduced iron group line blanketing) 
may provide signatures that the stellar envelope was in- 
deed of very low metallicity. 

The recently publicized SN 2007bi has been sug- 
gested to be the first convincing detection of a PI SN. 
iGal-Yam et al.l (|2009f ) compared the observed light curve 
of SN 2007bi to the models presented here, and found 
very good agreement with a helium core explosion of 
mass 100-110 Mq. We have shown here that the spec- 
trum of the helium core model is also in reasonable agree- 
ment with the observations, and in particular with the 
presence of strong magnesium lines. SN 2007bi is there- 
fore a compelling candidate for a PI SN, however other 
possib le scenarios have been suggested. iMoriva et al.l 
()2010[) have shown that the light curve could be ex- 
plained as well by the core collapse of a massive 43 Mq 
C/0 core. Another possibility is that the light curve was 
powered not by ^^Ni decay, but by energy injection from 
a highly magnetized remnant neutron st ar (a magnet ar 
iKasen and BildstenI IMot I WooslevI l20Tol ) ) . The distin- 
guishing feature of the pair SNe model is its exception- 
ally long rise time, but unfortunately the observations of 
SN 2007bi lacked data before the peak . At present, one 
appealing feature of the pair SNe model for SN 2007bi is 




Fig. 17. — Detectability of the shock breakout transient from pair instability supernovae as a function of redshift. The figure plots 
the observed AB magnitude (for different wavelength bands) of three red-supergiant models at the peak of the bolometric breakout light 
curve. (Note an AB magnitude of 28.9 corresponds to 10 nj.) Solid lines denote the redshift range for which the rest frame flux is redder 
than 1216 A (the Lyman alpha line); at higher redshifts, the flux is likely absorbed by intervening neutral hydrogen (dashed lines). The 
breakout transients of the blue supergiant models will be significantly dimmer than the models shown here, as the luminosity is lower and 
the spectrum peaks at much lower wavelengths. 



that it correctly predicts the explosion energy, given the 
progenitor star mass. In the other models, the energy is 
a free parameter that must be input by hand to match 
the observed light curves and spectra. 

Other luminou s observed sup ernovae, such 
as SN 2006gy (ISmith et all 12009) . SN 2005ap 
jpuimby e£an |2007D, and SN 2008es (|Gezari et al.l 
[2009) have occasionally been suggested to be examples 
of PI SNe. However, the light curve duration of these 
events (~ 50 — 100 days) are simply too short to be 
explained by the classic ^^Ni powered pair explosions 
explored here, which all last of order ~ 300 days. The 
pulsational pair scenario, which occurs for stellar masses 
just below those considered here, provide a possible 
explanation for these observations, as does magnetar 
energy injection. 

It is often lamented that pair SNe cannot be discovered 
as transients at high redshift - their intrinsically slow 
1 year) light curve evolution, prolonged by a (1 + 
z) cosmological time dilation, will exceed the lifetime of 
most observational surveys. Actually, the situation is 
not all that bad. For example, in our 100 Mq helium 
model the bolometric light curve declines (after peak) at 
at a rate of about 0.01 mag per day. At a redshift of 
z = 10, this results in a ~ 0.3 mag variation over an 
observer frame year. That may be within the sensitivity 
of future surveys with multi-year baselines. Moreover, we 
find that the decline is more rapid in the bluer bands (due 
to the increasing onset of iron group line blanketing at 



these wavelengths) such that variations in the restframe 
U- and B-bands are a factor of 2 — 3 larger. Even if the 
temporal variations are too hard to measure, the colors 
of PI SNe might be distinct enough that one could use 
them to select candidates for further follow up. 

An alternative approach for finding PI SNe would be 
to look for the brief, but very luminous emission of 
shock breakout. For the larger red-supergiant explo- 
sions, the breakout transients are 20 times as luminous 
(Lpcak ~ 10'*^ ergs s~^) as those of typical Type IIP su- 
pernovae, and last about four times as long ('^ 2 hours). 
The emission peaks in the rest-frame far ultraviolet (A ^ 
80- 170 A). Sensitive surveys in the infrared (2 — 5 /im) 
with a high cadence and a wide-field of view could, in 
principle, detect these transients out to high redshifts. 
Once detected, the subsequent ^^Ni powered light curve 
and spectra could be monitored, with some leisure, for 
many years to come. 
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